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hnd that the model with Type-II Yukawa interactions can better explain the experimental 
data among all. In this scenario, the couplings of the light CP-even Higgs boson h with 
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of tan p. In particular, we scrutinize a well-constrained region previously ignored by other 
analyses and study the phenomenology of the extra Higgs bosons at the Large Hadron 
Collider. 
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I. INTRODUCTION 

In the standard model (SM) of particle physics, the Higgs field plays a central role in breaking 
the electroweak symmetry and giving masses to other elementary particles. A direct consequence 
is the existence of a spin-0 Higgs boson that interacts with the other particles with strengths 
proportional to their respective masses. A new particle, denoted by h, with mass about 125.5 GeV 
has been recently discovered at the Large Hadron Collider (LHC) by both the ATLAS and CMS 
Collaborations [1|, 12|> an d is found to be compatible with the SM Higgs boson in the production 
rates of various channels. Moreover, the diphoton decay of h suggests that its spin is an even 
integer Q], and a preliminary angular analysis of the h — > ZZ* — > Al decay data prefers that the 
new boson has zero spin and positive parity j^]. 

Even though the new particle h is generally consistent with the SM Higgs boson, one obvious 
question is whether it is the only spin-0 particle in the SM or actually belongs to a larger scalar 
family. There are many well- motivated models with an extended Higgs sector, among which the 
two-Higgs doublet model (THDM) is one of the most popular and extensively analyzed classes. 
Such a structure of two Higgs doublet fields is also required for the Minimal Supersymmetric 
Standard Model (MSSM). By introducing a cousin doublet of the SM Higgs doublet, the model 
preserves the custodial symmetry in the kinetic term to keep the electroweak p parameter at unity 
at tree level. 

Recently, there are many analyses about the THDM's j^], with particular emphasis on the 
study of enhancement in the h — > 77 mode as the data show. A survey of generic models with 
an extended Higgs sector, including the THDM's, has been done by the authors jfl to show the 
correlation between the 77 and Zj modes. Obviously, a global analysis to the available data is 
required to disentangle among the possible new physics candidates 7|. Various scenarios in the 
THDM have been analyzed by using the Higgs boson search data at the LHC before the Moriond 



conference 



and those afterwards Q- 



In the THDM's, there may be tree-level flavor-changing neutral current (FCNC) interactions 
due to the mediations of neutral scalar bosons as both doublet fields can generally couple to the up- 
type and down-type quarks and charged leptons. There are several ways to avoid such dangerous 
FCNC processes, e^o., imposing a softly broken Z2 symmetry or assuming alignments in the 
Yukawa matrices 191 ] . In the former case, there are four independent types of Yukawa interactions, 



1 The model with an exact Z2 symmetry is known as the inert doublet model [li 
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depending on the charge assignments of the SM fermions under the Zi symmetry [20J. They are 
dubbed the Type-I, Type-II, Type-X and Type-Y THDM's 



22] . After electroweak symmetry 



breaking, all these scenarios contain the same physical Higgs bosons under the assumption of CP 
conservation in the Higgs sector: two charged Higgs bosons, two CP-even Higgs bosons, and one 
CP-odd Higgs boson. 

The Type-II THDM can best fit the latest Higgs search data from the LHC, with the preferred 
values tan/3 ~ 10 where tan/3 is defined as the ratio of the vacuum expectation values (VEV's) 
of the two doublet fields and SM-like couplings among the lighter CP-even Higgs boson and weak 
gauge bosons gnw 

However, we find that there is another parameter region that cannot be 
ruled out by the data at even 68% confidence level (CL) and is still consistent with the constraints 
of vacuum stability and perturbative unitarity. This region corresponds to the case where the 
values of the bottom and tau Yukawa couplings are smaller and have an opposite sign to their SM 
values. Thus, the Higgs to diphoton decay branching fraction is enhanced in comparison to the 
SM prediction, so that there is a possibility to explain the currently measured signal strength of 
the diphoton mode with a central value larger than 1.6 at both the ATLAS and CMS. Moreover, 
this region, ignored by previous analyses, has larger deviations (> 0(1%)) in the couplings guw 
from the corresponding SM values. This in turn has effects on the productions and decays of the 
other Higgs bosons at colliders. In this paper, we study the Higgs phenomenology in this part of 
the parameter space at the LHC. 

This paper is organized as follows. In Section HH we review the THDM with a Zi symmetry 
and classify four types of Higgs interactions with SM fermions. A soft Z<i symmetry breaking 
term is also introduced. In Section IIIIl we perform \ 2 fits of the four scenarios to the current 
data, singling out the Type-II interactions as the preferred one. We then find two separate regions 
in the model parameter space, noting that one of them has been ignored in previous analyses. 
We then concentrate on the region where the gnw couplings have larger deviations from the SM 
expectations, and study the collider phenomenology of the heavier Higgs bosons in Section IIV1 
After working out the branching fractions of both heavy CP-even Higgs boson, H, and CP-odd 
Higgs boson, A, we examine their single and pair productions at the LHC and compare the results 
with current search limits. Our findings are summarized in Section IVl 
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TABLE I: Charge assignments of the Z2 symmetry and the corresponding £/ factors in different scenarios 
of the THDM. 

II. THE TWO-HIGGS DOUBLET MODEL 



The THDM contains two SU(2) L doublet Higgs fields $1 and $2 with hypercharge Y = +1/2. 
In general, both doublet fields can couple to the SM fermions at the same time and induce FCNC's 
via the mediation of a scalar boson at tree level. To avoid such FCNC's, a softly broken discrete 
Z2 symmetry is introduced to the model, under which the doublet fields transform as $1 — > +$>i 



and $5 



-<I>2- The two doublet fields can be parameterized in the so-called Higgs basis as 




R{P) 




, with R(0) 




(1) 



where 



4> 



w 



4> 



and tan/3 = <$°)/<$°), 



(2) 



H+ 

with the VEV v = \J + (^Sj) 2 = 246 GeV. In the above parameterization of the component 
scalar fields, and z are the Nambu-Goldstone bosons becoming the longitudinal components of 
the and Z bosons, respectively. The physical charged Higgs boson (the CP-odd Higgs boson) 
is denoted by {A) and the CP-even Higgs bosons are expressed as hy and h' 2 . In general, (h^, 
h' 2 ) are not mass eigenstates and can mix with each other. The mass eigenstates of the CP-even 
Higgs bosons can be defined by introducing the mixing angle a as 2 



h' 



R(a-P) 




(3) 



We will assume that h is the newly discovered particle and is the lighter mass eigenstate; i.e., h is 
considered as the SM-like Higgs boson, and H is the heavier one. 



2 The CP-even scalar states in the basis of (/ii,/i2) can be directly related to those in the basis of (H,h) by 
(h u h 2 ) T = R(a)(H, h) T . 
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In the Higgs basis, the Yukawa Lagrangian with the Z 2 symmetry are given as 
V2 



Q L M d ($ + ) d R + Q L M U $ + ^) u R + L L M e ($ + £ e ^) e fl + h.c 



(4) 



where (/ = u, d or e) can be determined when we specify the Z 2 charges of the quarks and 
leptons, and Mf is the diagonalized fermion mass matrix. The charge conjugation of the scalar 
fields are denoted as $ = ir 2 <I>* and ^ = it 2 \I'* , where r 2 is the second Pauli matrix. There are 
four independent charge assignments of the Z2 parity as summarized in TABLE HJ The Yukawa 
interactions are expressed in terms of the CP-even Higgs mass eigenstates as 



Cy = ~ E — { [sin(/3- a) +£/cos(/3 -«)]///> 



f=u,d,e 



+ 



[ cos(/3 - a) - sin(/3 - «)]//# + zSign(/)^/ 75 /^} 
^ Vud u{m d ^ d P R - m u ^ u P L ) dH + + ^^^-VP R eH + + h.c. 



(5) 



where i^ jj are the projection operators for left- and right-handed fermions, respectively, and 
Sign(/) = +1 (-1) for / = d and e (/ = u). 
We consider the Higgs potential 

V = m 2 |$i| 2 + m 2 |$ 2 | 2 ~ m|(${$2 + h.c.) 

+ ^Axl^xl 4 + ^A 2 |$ 2 | 4 + A 3 |$i| 2 |$ 2 | 2 + A 4 |$i$2| 2 + ^A 5 [(*t* 2 ) 2 + h.c], (6) 

where m 2 and A5 are generally complex parameters. However, we assume CP-conserving THDM 
and, therefore, m| and A5 are taken to be real. Two of the parameters in the above potential m 2 
and m 2 can be related to other parameters using the tadpole conditions: 



m 1 =m 3 tan/3 — — (Aicos /3 + Asin /J), 
v 2 

m 2 . = m 2 cot /3 — 2~(^ 2 sm2 + A cos 2 /3), 



(7) 
(8) 



where A = A3 + A4 + A5 . The masses of and A are 



= M 2 -y(A 4 + A 5 ), m\ = M 2 - v 2 X 5 , 



(9) 



where 



sin /3 cos /3 



(10) 



The masses of the CP-even Higgs bosons are calculated as 

m \ = cos 2 (a - /?)M 1 2 1 + sin 2 (a - /5)M| 2 + sin 2(a - /3)M 2 2 , (11) 

m 2 h = sin 2 (a - /3)M 1 2 1 + cos 2 (a - /3)M 2 2 2 - sin 2(a - /3)M 2 2 , (12) 

where Af 2 l5 M 22 and M 2 2 are the elements of the mass matrix in the basis of (h^, h' 2 ) expressed by 

v 2 - 

M\ x = v 2 {\ x cos 4 (3 + A 2 sin 4 /3) + y Asin 2 2/3, (13a) 

M 2 2 2 = M 2 + v 2 sin 2 cos 2 /3(Ai + A 2 - 2A), (13b) 

2 2 

M\ 2 = V — sin 2/3(-Ai cos 2 /3 + A 2 sin 2 /3) + y sin 2/3 cos 2/3A, (13c) 



and the mixing angle is given as 



2M 2 



We note that the decoupling limit can be taken by letting M 2 — > oo 23|]. Among M 4 2 only M 2 2 
depends on M 2 , so that Af| 2 goes to infinity in this limit, corresponding to sin(a — j3) — > —1. 

It may be useful to write down the explicit formula for the hH + H~ vertex, which is important 
when we consider the loop contribution to the h — > 77 decay, as 

\ hH+H - = — 1— [2cos(a + /3)(m 2 - M 2 ) + sin(/3 - a) sm2(3(2m 2 H+ - m 2 )] . (15) 
u sin 

We note that when the sign of X hH + H - is negative (positive), the loop contribution to h — > 77 
has constructive (destructive) interference with the W boson loop contribution. 
The kinetic terms of the Higgs fields in the Higgs basis is given by 

£kin = |^$| 2 + l^l 2 - (16) 

where the covariant derivative 

D, = d,- l -gT.A^-^g'B^ (17) 

From this equation, it is seen that the first term is identical to the kinetic term of the SM Higgs 
boson. Therefore, the weak gauge boson masses are derived in the same way as in the SM. The 
couplings among the CP-even Higgs bosons and the weak gauge bosons, however, can be different 
from those in the SM, gfyy, because of the mixing between the two CP-even states. Explicitly, 

9hV V = 9hVV si 11 !/ 3 - a ) » 9HVV = 9hVV cos (/3 ~ a ) ■ ( 18 ) 
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Mode 


Data (ATLAS) 


/ £: 7 TcV+8 TeV (ATLAS) 


Data (CMS) 


JC: 7 TcV+8 TcV (CMS) 


h —> 77 


1 fi c+0.24+0.25 r 24 i 
l.OO_ 24 _ 18 


4.8 ftr 1 + 20.7 fb" 1 


0.77 ±0.27 (MVA) [28] 


5.1 fb" 1 + 19.6 fb" 1 








1.11 ±0.31 (Cut Based) [29] 


5.1 fb" 1 + 19.6 fb" 1 


h^fZZ 


1.5 ±0.4 [25] 


4.8 ftr 1 + 20.7 fb" 1 


0.91±£; 3 24 [30] 


5.1 fb- 1 + 19.6 fb- 1 


h -> WW 


1.01 ±0.31 [26] 


4.8 ftr 1 + 20.7 ftr 1 


0.76 ±0.21 [31] 


4.9 fb" 1 + 19.5 fb" 1 


h^bb 


-0.4 ± 1.0 [27] 


4.7 ftr 1 + 13 ftr 1 


i.3±8:S [32] 


5.0 fb" 1 + 12.1 fb" 1 


h — ^ TT 


0.8 ±0.7 [27] 


4.6 fb" 1 + 13 fb" 1 


1.1 ±0.4 [33] 


4.9 fb" 1 + 19.4 fb- 1 



TABLE II: Signal strengths of the five processes measured at the ATLAS and CMS. 



III. DATA FITTING 

So far, the Higgs boson search at the LHC has been done in the following five processes; pp —> 77, 
pp ZZ*, pp — > WW*, pp — > tt, qq' — > Vh — > Vbb, where pp — > h and qq' — > Vh indicate the 
inclusive Higgs boson production and the vector boson associate production, respectively. The 
signal strength for each of the channels is defined as 

Rcf _ <tf* X BR(h -> X)™ 

^ X ~ xBR(/w!)SM' 1 ' 

where of ef [o^ M ] and BR(h -> X) Rc{ [BR(h -> X) SM ] are the reference value [SM prediction] of the 
Higgs production cross section and that of the branching fraction of the h — > X decay, respectively. 
Experimental data of fix are listed in TABLE |TTJ The average signal strength of the ATLAS and 
CMS data can be calculated using 

cxp = /£ ATLAS X / ,ATLAS + /£ CMS x/i CMS 
— J £ ATLAS _|_ J £CMS ' 

where J £atlas (/ £cms) is the integrated luminosity of the ATLAS (CMS) Collaboration. We 
then obtain the average signal strengths as 

^ = 1.22 ±0.31, /4 X | = 1.21 ± 0.35, ^ w = 0.89 ± 0.27, 

^ = 0.44 ± 0.87, {if? = 0.97 ± 0.55, (21) 

where, to be conservative, the standard deviation of each signal strength is derived by using the 
larger error when the error bar is asymmetric. For /x^* p , we used the experimental value based on 
the MVA method. With the input of Eq. (12 1 H . a x 2 value can be calculated for each reference value 
as 

/„ cx p _ „Rcf\ 2 
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Model 


v 2 . 

A. mm 


tan j3 


SM 


1.45 


- 


Type-I 


1.47 


3.1 


Typc-II 


1.21 


4.3 


Type-X 


1.45 


2.3 


Type-Y 


1.24 


4.5 



TABLE III: Minimal x 2 value (x^ in ) for the SM and the THDM's with m$ = M = 300 GeV and 5 = 10~ 4 . 



where A^t^ p is the one standard deviation of fJ°x P - 

We evaluate the x 2 values for the four scenarios (Type-I, Type-II, Type-X and Type-Y) defined 
in the previous section. In addition to the VEV v = 246 GeV and the mass of the SM-like Higgs 
boson mh = 126 GeV, there are six free parameters in the Higgs sector: the three masses for 
the extra Higgs bosons mjj, t>^a and m H +, the two mixing angles a and /3, and the dimensionful 
parameter M 2 representing the scale of softly broken Z2 symmetry. In the following discussion, 
we consider the case where the masses of the extra Higgs bosons are taken to be degenerate, 
m$ = mjj = niA = m H+i f° r simplicity and to avoid additional contributions to the electroweak 
p parameter. In addition a large difference between m$ and M means a large value of Higgs 
qurtec couplings, it can be excluded by the constraints from the perturbative unitarity and also 
the vacuum stability, so that we take M to be the same value as m$. The top Yukawa coupling 
is constrained by perturbativity to be \yt\ 2 < 4-7T with yt = V^rnt / (v sin /3) . This gives the lower 
bound of tan /3 > 0.3. 

We here summarize constraints from B physics studies. First, the B — > X s ~f data demand that 
the mass of the charged Higgs boson, mjj+ , to be less than 295 GeV at 95% CL 34| in the Type-II 



and Type-Y THDM's when tan/3 > 2. When tan/3 < 2, m H + < 300 GeV is excluded at 
CL in all types of THDM's 35J. A similar but slightly weaker bound for tan/3 with m H + < 300 



GeV has also been given by the data of the Z — >■ bb decay 



351 ]. From the B u — > tv, B — > Dtu 



and K — > ny data, the Type-II THDM with tan /3 > 30 is excluded at 95% CL for mjj+ = 300 



GeV 



35|, 



36J. In accord with the above-mentioned constraints, we will take m<j> = 300 GeV and 



tan (5 > 2 in the following numerical analyses. 

We introduce a parameter 5 = 1 — sin(/3 — a) to describe the deviation in the hVV couplings 
from the corresponding SM values. Therefore, the mixing angle a is determined for a given pair of 
5 and tan/3. The deviations in the Higgs boson couplings with the weak gauge bosons c^vv an d 
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FIG. 1: /i™ DM as a function of tan/3 for the Type-I (upper left), Type-II (upper right), Type-X (lower left) 
and Type-Y (lower-right) THDM, with S = 1(T 2 and m$ = M = 300 GeV. 



in the Yukawa couplings Chff are 

c hVV = 1-5, for V = W and Z, (23) 
c w = l + C/\/27 + 0(J), for / = t, b and r, (24) 

where the factor £y (£ f = ^ = ^ and £ r = £ e ) is listed in TABLE [Ij 

In TABLE IIIIl the minimal value of x 2 (x m in) ano - * ne corresponding tan/3 value are listed for 
the SM and the THDM's with 5 = 1CT 4 and m$ = M = 300 GeV. The smallest value of x min 
among these models is obtained in the Type-II THDM with tan /3 = 10.2. Except for the Type-I 
THDM, the value of Xmin s ^ a Y s almost the same for 5 < 10~ 4 . But the value of tan/3 at the x 2 
minimum gets larger than that for 5 = 10~ 4 . In the Type-I THDM, the value of x m in approaches 
3.26, still larger than the SM value, in the limit of 8 — > with tan /3 ~ 4.3. This is due to a 
destructive loop contribution to the decay rate of h — > 77 with the W boson loop contribution 
(see Eq. ([Pol)). 
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FIG. 2: ^™ DM as a function of tan/? for the Type-II THDM with m # = M = 300 GcV. The upper-left, 
upper-right, lower-left, and lower-right panels show the results for 5 = 0.1, 10~ 2 , 10~ 3 and 10~ 4 , respectively. 



In Fig. [lj values for //™ DM are plotted as a function of tan /3 in the Type-I (upper- left), Type-II 
(upper-right), Type-X (lower-left) and Type-Y (lower-right) THDM's with 5 = 10~ 2 . It is seen that 
in the Type-II THDM, both /^™ DM and /i™ DM are around 1.5 while all the others approach in 
the large tan (5 region. This is because both bottom and tau Yukawa couplings are enhanced by the 
same factor 1 — tan Py25, while the top Yukawa and gauge couplings are almost the same as their 
SM values. Thus, the Higgs production cross section is barely changed from the SM prediction. 
Yet the branching fractions of h — > rr and h —> bb modes are as large as about 10% and 90%, 
respectively, roughly 1.5 times larger than the SM values. On the other hand, in the Type-X (Type- 
Y) THDM, only the tau (bottom) Yukawa coupling can be enhanced by the factor of 1 — tan (3^25, 
so that only /i™ DM (/^™ DM ) is larger than 1 in the large tan/3 region. The asymptotic values 

TH 
l bb 

tan/3 limit. 



of /i™ DM (^™ DM ) are about 14 (1.7) in the Type-X (Type-Y) THDM, respectively, in the large 
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We note that all the //£ HDM values are around 1 at tan/3 ~ 15 where c^bb and c^ TT (Type-II), 
Chrr (Type-X) and Chbb (Type-Y) have an opposite sign to their corresponding SM couplings, while 
the other cuff's and Chvv are almost the same as the SM values. In the Type-I THDM, all the 
Yukawa couplings are modified by the same factor, and they are close to the SM values in the 
large tan/3 region so that all the /U™ DM values approach to 1 except for /i™ DM . Only the value 
of /i™ DM becomes smaller than 1 for large tan/3 because of the destructive contribution of the 
charged Higgs boson loop 3 . 

In Fig. [2j values of /U™ DM are displayed as a function of tan (3 in the Type-II THDM with 
5 = 0.1 (upper-left), 10~ 2 (upper-right), 10 -3 (lower-left), and 10~ 4 (lower-left). It is seen that 
in the cases with smaller values of 5, the value of /i™ DM (At™ DM ) gets quite close to that of 
^™ DM (/i™ DM ), because the decay rates of h — > bb and h — > tt can be changed by the same 
factor c\ bb = (? hrr = (1 — tan/3\/25) 2 , and the other decay rates (h — > 77 and h — > VV) and the 
production cross sections (the gluon fusion and the gauge boson associated production) are almost 
the same as in the SM. We observe a peak for both /i™ DM and /i™ DM , especially for the case 
with 5 < 0.1. This is caused by vanishing bottom and tau Yukawa couplings. The location of the 
peak can be determined by solving Chbb = Qm-t = 0, from which tan (3 ~ 1 /\^25. 

In Fig. [3 x 2 is plotted as a function of tan /3 for 5 = 0.1, 10^ 2 , 10^ 3 and 10~ 4 in the THDM's of 
Type-I (upper-left), Type-II (upper-right), Type-X (lower- left) and Type-Y (lower-right). In the 
Type-II and Type-Y THDM, there is a valley in the curve for each case of 5 except for 5 = 10~ 4 , 
where the bottom of the valley gives the minimal x 2 value. The value of tan /3 corresponding to the 
bottom of the valley is slightly smaller than the one for all the /^™ DM values to be unity, as seen 
in Fig. [TJ In the Type-X THDM, there is a similar valley in each value of S except for 5 = 10~ 4 , 
yet the bottom of that does not correspond to the global x 2 minimum. This is because the h — >■ 77 
channel is not sufficiently enhanced in comparison with that in the Type-II and Type-Y THDM's 
(see the lower-left plot in Fig. [1]). In the Type-I THDM, there is no such a valley, and the minimum 
of x 2 is obtained at around tan/3 = 3.1 independent of the value of 5. 

Fig. [J] shows the contour plot for x 2 m the tan/3-<5 plane in the Type-I (upper-left), Type-II 
(upper-right), Type-X (lower- left) and Type-Y (lower-right) THDM's. The cross X indicates the 
point that gives the minimal x 2 value, Xmin- ^ ne rec ^' oran g e an d black curves are respectively the 
contours corresponding to Xmi n + 0.2, Xmi n + 0.5 and Xmin + 2.3, the last case corresponding to 68% 



3 The same effect occurs in the other three modeis as weii. But it is masked by the iarge branching fractions of 
h — > ff modes; e.g., large branching fractions of h — > bb and h — > tt in the Type-II THDM. 
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Type I Type II 




tan(3 tan(3 

FIG. 3: x 2 as a function of tan/3 for several values of 5 and m$ = M = 300 GcV. The upper-left, upper- 
right, lower-left and lower-right panels show the results in the Type-I, Type-II, Type-X and Type-Y THDM, 
respectively. 



37| 



CL. The light blue (dark magenta) regions are excluded by the constraint of vacuum stability 
(perturbative unitarity [38j]). In the upper-right figure, the light green area shows the excluded 
region by the B physics studies. It is seen that there is an isolated narrow region consistent with 
data at 68% CL for Type-II, Type-X and Type-Y THDM's, and is not excluded by the vacuum 
stability and perturbative unitarity conditions. The values of tan (3 and 5 in this region correspond 
roughly to the bottom of the valleys in Fig. [3j 

From the above analysis, we here conclude that Type-II THDM can best explain the current 
LHC data among all the THDM's with softly broken Z2 symmetry. This is because the branching 
fraction of the h — > 77 mode can be enhanced due to the suppressed decay rates of h — > bb and 
h — > tt modes. There are two regions where the data can be well reproduced in the model of 
Type-II. One of them is the region with a very small 5; e.g., 5 < 0.01. The other is the region 
satisfying Chbb — Chrr — — 1 j which gives the relation 5 tan /3 ~ y/25. This is a good approximation 
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FIG. 4: Contour plots of x 2 values in the tan /3-d plane with to$ = M = 300 GeV. The upper-left, upper- 
right, lower-left and lower-right panels show the results for the models of Type-I, Type-II, Type-X and 
Type-Y, respectively. The coordinate marked by X has the minimum \ 2 value. The red, orange and black 
curves are respectively contours for xLin + 0-2, Xmin + 0-5 and xLin + 2-3 (68% CL). The light blue (dark 
magenta) regions are excluded by the constraint of vacuum stability (perturbative unitarity). The light 
green region in the upper-right panel is excluded by the B physics studies mentioned in the main text. 



especially for small values of 5. 

IV. PHENOMENOLOGY OF EXTRA HIGGS BOSONS 

In this section, we discuss the phenomenology of the extra Higgs bosons, i.e., the heavier CP- 
even Higgs boson H, the CP-odd Higgs boson A and the charged Higgs bosons in the parameter 
regions favored by the LHC data. In the previous section, we found that the LHC data could be 
better explained in the Type-II THDM compared to the other THDM's. Therefore, we focus on 
this scenario with m$ = M = 300 GeV as in the previous section. Furthermore, when 5 is larger 
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FIG. 5: Decay branching fractions of H as a function of 8 (left panel) and tan/3 (right panel) for to$ = 
M = 300 GcV in the Type-II THDM. The band for each curve reflects parameter uncertainty at the 68% 
CL as determined by the Higgs data. The light-colored parts are excluded by the LHC search data [42 [. 



than 10~ 2 , the favored parameter space on the tan @-8 plane is restricted to a narrow band, so 
that the value of tan/3 can be approximately determined for each given value of 5. Such larger 
deviations in the hZZ and hWW couplings (5 > 10 -2 ) are expected to be measurable at future 



39]. 



MUM for the MSSM whose Higgs 



linear colliders such as the International Linear Collider 

The extra Higgs boson searches have been reported in Refs. 
sector corresponds to that of the Type-II THDM with a super symmetric relation among the Higgs 
parameters. So far, the search for the charged Higgs bosons has only be done in the — > t^v and 
— > cs channels via the top quark decay t — > H + b. When the charged Higgs boson mass is larger 
than the top quark mass, no bound on its mass is currently available at the LHC 4CJ, |4l|]. On the 



other hand, the search for the extra neutral Higgs bosons H and A has been done in the gluon fusion 



fig 



bbH/A -»• bbrr. 



4J by the CMS and in 



process: gg — > H/A — > tt/h/i and the bottom quark associate processes: 
Current excluded regions in the m^-tan (3 plane have been shown in Ref. 

Ref. Q h y the ATLAS. F„ r eWe. if mA is taken to be 300 GeV, the u ppe r hound tor the 
value of tan/3 is given by 7.58 |42l |. We note that this bound for tan/3 can be modified in the 
(non-supersymmetric) Type-II THDM, for the production rates for H and A and also the decay 
branching ratios of H/A — > tt can be different from those in the MSSM. 

First, we evaluate the decay branching fractions of H, A and H^. Fig. [5] shows the decay 
branching fractions of H as a function of 5 (left panel) and tan/3 (right panel) in the Type-II 
THDM. The band for each curve indicates the uncertainty in the parameter choice at the 68% 
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FIG. 6: Same as Fig.0 but for A. 
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FIG. 7: Same as Fig. El but for d 



CL. Note that we have restricted ourselves to the region with 5 > 1(T 2 and 2 < tan (3 < 10 in the 
upper-right panel of Fig. HI T he parts with light colors are excluded by the heavy neutral Higgs 
boson search at the LHC 42|. It is seen that the gauge boson pair decays of H (H — > WW and 
H ZZ) become more important than the fermionic decays (H — > bb and H — > tt) when the 
value of 5 (tan/3) is taken to be larger (smaller). The situation where H mainly decays into the 
gauge boson pairs does not happen in the MSSM, for 5 is suppressed in the large tua regime. For 
example, when tua is taken to be 300 GeV, the value of 5 is smaller than 10~ 2 [44], so that only 
the fermionic decays of H dominate. 

Fig. [6] shows the decay branching fractions of A as a function of 5 (left panel) and tan (3 (right 
panel) in the same setup as in Fig. [SJ The decay rate of A — > hZ (A — > bb) is proportional to 
cos(/3 — a) 2 ~ 25 (tan 2 /3), so that the magnitude of the branching fraction of A — > hZ (A — > bb) 
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FIG. 8: Gluon fusion production cross sections for H and A in units of fb as a function of 6 (left panel) and 
tan j3 (right panel) for to$ = M = 300 GcV in the Type-II THDM. The collision energy is assumed to be 
8 TeV. The band for each curve reflects parameter uncertainty at the 68% CL as determined by the Higgs 
data. The light-colored parts are excluded by the LHC search data 42 1 . 



increases (decreases) when 5 gets larger values. When 5 > 0.03 or tan (3 < 7, the branching fraction 
of A —7- hZ is larger than that of the A — > bb decay. 

In Fig. [7J the decay branching fractions of are shown as a function of 5 (left panel) and 
tan/3 (right panel) in the same setup as in Fig. In the region of 5 > 10~ 2 , the H + — > tb and 
— > hW^ modes are dominant. 

Next, we discuss the production of the extra Higgs bosons. The extra neutral Higgs bosons H 
and A are mainly produced via the gluon fusion process: gg — > H/A. The production cross section 
is given by 

a(gg -> H/A) = a(gg -> h SM ) x ^ % (25) 

m H/A T ( h SM -> 99) 

where cr(gg — > /ism) is the gluon fusion production cross section of the SM Higgs boson hsyi and 
T(/isM — > 99) is the decay rate of hsM — > 99- We note that the vector boson fusion production 
mechanism is not useful to produce H and A because the cross section for H is proportional to 
cos(/3 — a) 2 ~ 25 and that for A is zero at the tree level due to the absence of the VVA vertex. 

In Fig. [51 the gluon fusion production cross sections for H and A are shown as a function of 5 
(left panel) and tan /3 (right panel) for the collision energy to be 8 TeV in the Type-II THDM. The 
cross sections are calculated using t he g luon fusion cross section of a fictitious 300-GeV Higgs boson 
in the SM, whose value is 3.606 pb 



451 ] . The band for each curve reflects parameter uncertainty at 



the 68% CL as determined 



3y the Higgs data. The predictions shown by light colors are excluded 



by the LHC search data [421 ]. Except for 5 ~ 10 2 , the cross section of A is larger than that of H 
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FIG. 9: Same as Fig. [51 but for the bottom quark associate production cross section. 



because the Hti coupling is proportional to ~ (v 2(5 — cot /3), so that larger values of 5 give smaller 
coupling constants, while the Att coupling is simply proportional to cot /3. No dependence of 5 in 
the Aff couplings can also be seen in the cross section of A in the right plot of Fig. [HJ namely, the 
gg — >• A cross section is given by a curve without a band. 

As another single production mechanism for H and A, the bottom quark associate processes: 
gg — > bbH/A can be important, especially in the case of large tan/3. The cross sections of these 
processes are proportional to (V25 + tan/3) 2 and tan 2 j3 for H and A, respectively, in the Type-II 
(and Type-Y) THDM. 

In Fig. [9J the bottom quark associate production cross sections for H and A are plotted as a 
function of 5 (left panel) and tan f3 (right panel) in the same setup as in Fig. The cross sections 
are calculated by scaling from that with tan (3 = 1 and m$ = 300 GeV whose value is obtained 
as 5.34 fb using the MadGraph package 



46| and the CTEQ6L parton distribution functions (PDF's). 
The cross section for H is almost the same as that for A as they have almost the same scaling 
behavior, and the maximum value is about 300 fb when tan /3 is taken to be 15. 

Regarding the charged Higgs boson production, the gb — > H^t process is important, especially 
when the charged Higgs boson mass is larger than the top quark mass. The tbH^ vertex is 
proportional to m^tan (3 + mt cot /3, so that the cross section reaches the minimum when tan /? = 
yjm t /m b ~ 7.5. 

In Fig. [TUl the gb — > H^t cross section is shown as a function of 5 (left panel) and tan /3 (right 
panel) in the same setup as in Fig. [HJ The cross section is calculated by scaling from that with 
tan j3 = 1 and m$ = 300 GeV whose value is obtained as 120 fb using the MadGraph package and 
the CTEQ6L PDF's. The cross section has a minimum of ~ 7 fb when tan (3 ~ 7.5. The maximum 
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FIG. 10: Same as Fig. [51 but for the gb H ± t process. 



BF B(H -> W+W 


") B(H->ZZ) 


B(A->hZ) 




6(i7 ± -> t&) 


[%} 65-68 


29-30 


95-97 


45-53 


47-55 


CS a(gg -> ff) 




cr(.gg — > bbH) 




a(gb -> if"*) 


[fb] 211-300 


1285-1921 


32-43 


29-42 


19-26 



TABLE IV: Branching fractions (BF) and cross sections (CS) of if, A and H ± for 5 = 0.2 and m$ = M =300 
GeV. The allowed range of tan j3 at the 68% CL, 2.3 < tan/3 < 2.8, is reflected in the range of each quantity. 

value of the cross section is about 33 fb when the tan f3 ~ 2, corresponding to 5 ~ 0.3. 

There are also pair production processes for the extra Higgs bosons via the s-channel gauge 
boson mediation such as qq -)• j*/Z* ~* H + H~, qq ->■ Z* -»■ .£L4, qq' ->■ TF ± * ->■ iliT* and 
qq' —7- VF 1 * 1 * —7- AH^ . However, the cross sections are suppressed with increasing masses of the 
extra Higgs bosons. For example, when m$ = 300 GeV, the cross section of pp — > AH + is about 
1.5 fb at the collision energy of 8 TeV, and the other cross sections are even smaller [44]. 



Finally, we would like to discuss the signal events for the extra Higgs bosons at the LHC in the 
favored parameter region. As seen in Figs. El [6] and [3 H (A) mainly decays into the gauge boson 
pairs W + W~ and ZZ {h and the Z boson), while mainly decay into hW^ and tb in the large 
5 region (5 > 10 -2 ). As an example, we here consider the case with 5 = 0.2 in which the allowed 
range of tan f3 at the 68% CL is 2.3 < tan f3 < 2.8. 

TABLE IIVI lists the branching fractions and production cross sections of H, A and . For 
H (A), the gluon fusion production cross section is about 7 (50) times larger than that of the 
corresponding bottom quark associate production process. The gb — > H^t production cross section 
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Events £+£-fi T i+i'l+l' bbl+i' £+£-bbbf] T 
Signal [fb] 7.1 - 9.5 0.29- 0.40 24- 37 0.29 - 0.34 
Backg. [fb] 1.73 x 10 3 22.5 67.3 0.048 



TABLE V: Cross sections of the signal and background processes for the collision energy of 8 TeV at the 
LHC, obtained using the MadGraph package and the CTEQ6L PDF's. 

is about 19 — 26 fb. The following final states of the signal processes are important to test the 
scenario with large values of 5: 

gg^ H W + W~/ZZ -> $ T /£ + r £ + T for H, (26) 
gg -> A -> hZ -> bU + £~ for A, (27) 
gb^ H~t^ hW~bW + ^ £ + rbbb$ T for H ± , (28) 

where ^ represent or /r^. The signal processes gg — >■ -ff — >■ W + W~ jZZ have the main 
backgrounds from — >• W + W~ /ZZ, and are the same as the SM Higgs boson with a fictitious mass 
of 300 GeV. For the signal event gg — > A —> hZ (gb — > tH~), the pp — > ZZ and qq — >■ — >■ h$y[Z 
{gb — > W + W~hsMb) processes can be main backgrounds. 

In TABLE |V1 the signal and background cross sections are listed assuming the LHC collision 
energy to be 8 TeV. We here assume that the b-jets can be identified with 100% efficiency. Except 
for the signal events from H (first two columns), the other two cross sections for the signal events 
are comparable or larger than those of the background. The cross sections for the signal events 
from H is smaller by about two order of the magnitude than that from the backgrounds. However, 
for the final states of £ + £~ IpT and £ + £~£ + £~ coming from the H signal processes, the transverse 
mass distribution (471 ] and the invariant mass distribution of the £ + £~£ + £~ system can be useful in 
increasing the signal-to-background ratio. The Jacobian peak (resonance peak) of the signal in the 
transverse mass distribution (the invariant mass distribution) can be used to determine the mass 
of H, while such a characteristic feature cannot be observed in the background events. 

We now briefly comment on the Higgs phenomenology in the other parameter region preferred by 
the LHC data, i.e., the region with 5 < 10" 2 in Fig. H In that case, the HVV, AhZ and H^hW^ 
couplings are much suppressed, and the decay branching fractions of H — > VV, A — > hZ and 
are negligible. Instead of such gauge boson associate decay modes, these heavy Higgs 
bosons mainly decay into fermion pairs whose pattern depends on the type of Yukawa interactions, 
the masses of the extra Higgs bosons, and tan j3. For example, H and A can decay dominantly 
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into bb (rr) for large tan (3 in the Type-I, Type-II and Type-Y (Type-X) THDM when the masses 
of both H and A are smaller than 2m^. On the other hand, the main decay mode of can be rv 



(cs and cb) for large tan /3 in the Type-I, Type-II and 
Higgs boson mass is smaller than the top quark mass 



^pe-X (Type-Y) THDM when the charged 



21). 



V. CONCLUSIONS 



In this work, we have considered the two-Higgs doublet models (THDM's) with softly broken 
Z2 symmetry introduced to avoid flavor-changing neutral currents at the tree level and the four 
independent scenarios (Type-I, Type-II, Type-X and Type-Y) differing in the Yukawa interactions. 
We have scanned the parameter regions that can explain the current Higgs boson search data 
recently reported by the ATLAS and CMS Collaborations. We have found that the Type-II THDM 
can best explain the data among all, with two separate parameter regions on the tan (3-5 plane at 
the 68% CL. One region has smaller 5 (< 10 -2 ), and the other implies the relations Chbbi Ch TT ~ — 1. 
The latter case, with tan /3 ~ and 5 > 10~ 2 , is ignored in previous analyses. Based on this 

finding, we have then studied the phenomenology of the extra Higgs bosons in this region. The 
heavy CP-even Higgs boson primarily decays into a pair of W and Z bosons, different from the 
corresponding particle in the minimal supersymmetric standard model. The CP-odd Higgs boson 
mainly decays to the SM-like Higgs and Z bosons, and the charged Higgs boson to the SM-like 
Higgs and W bosons or a pair of top and bottom quarks. We have also computed the cross sections 
of various signal and background events for the production of the extra Higgs bosons. 
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Note added: 



481 ] in which a 



While this manuscript was being written up, we noticed the appearance of Ref. 
global fit to the Higgs boson search data at the LHC was discussed for various THDM's, including 
the softly broken Z<i symmetric THDM's considered in this work. 
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